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Abstract 

We propose a semianalytical ion dynamics model for a collisionless radio frequency biased 
sheath. The model uses bulk plasma conditions and electrode boundary condition to 
predict ion impact energy distribution and electrical properties of the sheath. The proposed 
model acco un ts for ion inertia and ion current modulation at bias frequencies that are of 
the same order of magnitude as the ion plasma frequency. A relaxation equation for ion 
current oscillations is derived which is coupled with a damped potential equation in order 
to model ion inertia effects. We find that inclusion of ion current modulation in the sheath 
model shows marked improvements in the predictions of sheath electrical properties and 
ion energy distribution function. 
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I. INTRODUCTION 


High density, low pressure plasma tools are widely used for etching of semiconductor 
materials. 1 These tools usually employ a radio frequency (rf) biased electrode on which 
the wafer is placed. This rf bias controls the impact energies of the ions arriving at the 
wafer, which can be tuned to attain desirable etch rate and selectivity. The ions are 
mostly produced by electron impact ionization in the center of the plasma. They are then 
accelerated toward the substrate and the reactor chamber walls by strong, time-varying 
electric field in the sub-millimeter scale sheath that is formed at the plasma boundary. 
Thus, it is important to compute the acceleration of the ions as they traverse the non- 
collisional sheath in order to predict ion energy distribution (IED) at the wafer. The 
dynamics of the ions in the sheath is governed by the inertial and electrical forces. The 
electrons, however, due to their low inertia and high thermalization are in equilibrium 
with the instantaneous potential and follow a Boltzmann distribution at bias frequencies 
of interest. 

In the absence of a time- varying component of the potential, the ion kinetic energy 
gain is simply balanced by the drop in the electrostatic potential in the sheath. However, 
when the rf bias is applied at the electrode, the sheath electric field oscillates in time. 
This makes a theoretical analysis more complex. Although, a numerical modeling of ion 
dynamics is possible, it is a formidable task to self-consistently simulate it as a part of 
a multidimensional reactor model. Hence, it is desirable to have a semianalytical model 
that allows us to determine the ion impact energies at the wafer and electrical properties 
of the sheath. This semianalytical model can then be used as a boundary condition in a 
full reactor simulation. 2 The expressions for sheath electrical properties are also useful in 
a plasma-circuit interaction models. 3 

A theoretical analysis depends on simplifications based on the relative values of the 
rf bias frequency, u>, and the ion plasma frequency, u p i (or the ratio (3 — u/upi), which 
determines how much the ions respond to the time variation of the sheath potential. If the 
applied rf bias frequency at the electrode is much different from the ion plasma frequency 
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the ion dynamics can be considerably simplified. In case of a small applied rf bias frequency 
compared to the ion plasma frequency (0 = u>/u P i <C 1), the ions have enough time to fully 
adjust to the changes in the sheath potential. This is the equilibrium condition where the 
ion density and velocity at a point in the sheath is only a function of the instantaneous 
value of the local potential. At the other extreme, when the applied rf bias frequency is 
much higher (0 = u/u pi » 1), the ions remain in a frozen condition as they are unable to 
adjust to the rapid changes of potential. Under this condition the ion velocity and density 
ran be determined by the local, time averaged value of the potential. Analytical models 
for the equilibrium and frozen conditions have been obtained by various researchers. 4 

In this paper we will focus on the intermediate frequency range when the rf bias fre- 
quency is of the same order of magnitude as the ion plasma frequency (0 ~ 1). In this 
nonequilibrium regime, the ions only partially respond to the time variation of the sheath 
electric field. This makes an analytical simplification of the ion fluid equations difficult. 
Riley 7,8 proposed a unified sheath model for rf driven sheaths that is applicable for all 
frequencies. In the intermediate frequencies, the ion inertia effects are modeled using a 
damped potential, V, which is obtained from a relaxation equation, 


dV V - V 

~dt ~ Ti ’ 


( 1 ) 


where t* is the ion transit time through the sheath. It was assumed that the ions follow 
the damped potential instead of the rapidly changing sheath potential, V . This equation 
is derived by spatially integrating the ion momentum equation. 7 In order to simplify the 
analysis and obtain a close form analytical expression for the electric field, the unified 
sheath model assumes a constant ion current through the sheath at all time. This assump- 
tion appears to be rather poor and simplifies the analysis at the cost of accuracy of the 
model. In Ref. 9, we showed that ion current through a sheath can vary by as much as 
70% and has significant effects on the predictions of the electric field, sheath capacitance, 
IEDs, etc. We then derived a correction term in the expression for the electric field which 
accounted for the spatio-temporal variation of ion current in the sheath. However, we used 
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an empirical approach in the above work to determine the magnitude of the ion current 
oscillations. This restricted the validity of the model to only cases when the sheath po- 
tential varies sinusoidally. In this work, we improve this model and propose a modified 
set of equations that does not depend on empirical expressions and is valid over a wide 
range of potential and current conditions. We propose a relaxation equation for ion current 
oscillations which must be integrated in time along with Eq (1). 

The rest of this article is organized as follows. In Sec. II, we present the continuum 
equations that govern ion transport in a non-collisional sheath and the relevant boundary 
conditions. In Sec. Ill, we describe the details of the theoretical analysis and the model 
formulation. In Sec IV, we discuss the results obtained and compare them with the unified 
sheath model results. Both models are evaluated against the direct numerical solution of 
the ion transport equations. 

II. GOVERNING EQUATIONS 

Ions in a non-collisional sheath can be assumed to be non-thermal with high directed 
kinetic energy. This allows us to use the cold ion fluid equations. 


dn dnu 

dt "** dx 
du du 


= 0 


e dV 


dt U dx m dx 

d 2y t 

77 = (n - n e ) 

dx 2 e 0 


( 2 ) 


n and n e are the ion and electron number densities, u is the ion velocity, and V is the 
potential in the sheath, e and m are the unit charge and ion mass respectively. The 
electrons sire assumed to follow a Boltzmann distribution, 


n e = no exp 



( 3 ) 


where T e is the electron temperature written in volts and no is the electron (and ion) 
density at the sheath-presheath edge, where a steady state condition is assumed. We 
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further assume that the electric field at this boundary is small compared to the fields 
in the sheath. This enables us to apply asymptotic boundary conditions as described by 
Riemann and coworkers. 10,11 For a numerical solution of Eq. (2), the location of the sheath- 
presheath boundary is immaterial as long as it is sufficiently far away from the electron 
front. The boundary condition at the sheath-presheath edge is based on quasineutrality 
and Bohm criterion. Therefore, at x = -L (sheath-presheath boundary), 

n = n e = n 0 , u = u B - and V = 0 (4) 

where, u B is the ion sonic velocity. A reference potential is chosen at this boundary. Since 
the ions attain supersonic velocity in the sheath, an extrapolation condition is used for n 
and nu at the electrode boundary, x = 0. In order to close the system, another condition 
is needed at the electrode boundary which is determined from circuit considerations. 

III. SEMIANALYTICAL MODEL 

A. Conservation Laws 

A particle or ion current conservation law can be obtained by integrating the ion conti- 
nuity equation [first equation in Eq. (2)]. It is evident that ion conduction current is not 
necessarily constant in space and time. 


enu — enou B = [ endx = J, (5) 

Ot J-L 

The above expression shows that the ion conduction current is dependent on the rate at 
which the accumulated ion charge in the sheath changes in time. This is expected to 
be limited by ion inertia. From Eq. (5), we note that ion current has two components; 
a constant ion current, eriQU B , entering the sheath from the plasma and a time varying 
component, is an unknown and will be determined later. 

The kinetic energy conservation of ions is obtained by integrating the ion momentum 
equations and defining a damped potential, V. 
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The derivation of the above equation is given in Ref. 7. The assumptions made in order 
to derive this equation are not affected when a time varying ion current is introduced in 
the analysis. The ion velocity is now conveniently written as 



The damped potential is determined by integrating Eq (1). 


(7) * 


B. Electric Field 

The above conservation laws can be used to write the Poisson’s equation in the following 
form 


&V 

dx 2 




2V\-i/2 


( 8 ) 


An expression for the electric field can be obtained by multiplying the Poisson’s equation 
by dV/dx and integrating over the sheath length. This yields, 


E 2 

2 


s 

2 


= h + h + 



( 9 ) 


where Eq is the electric field at the sheath-presheath edge and I\ and / 2 are the following 
integrals: 


h = 


h = ~ 


e 

en 0 


/>€)- 
if (=*=)(- 


dV\ 

2V'\~V 2 


( 10 ) 


where V' is a dummy variable in the above integral. We note that Eq (9) is the same 
expression that is obtained by Riley 7 with an additional correction term, / 2 . This term 
accounts for the ion flux oscillations. As expected, h vanishes when the ion current is 
assumed constant (t.e. J x = 0). 
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In order to evaluate the integrals, h and / 2 , certain simplifying assumptions regarding the 
spatial variation of V and are necessary. For V we adopt a linear relation, V{x, t ) = 
a(t)V(x,t), as proposed by Riley. 7 . a(t) is a function of time. For Ji, we assume the 
following relation, 


Mx,t) =7(t)[v(x,t)r <q (u) 

where 7 is a constant of proportionality and r is an exponent to be determined later. These 
quantities are only functions of time. The above relation is necessary in order to evaluate 
J 2 , its validity, however, can only be justified by the improved performance of the proposed 
model. The integrals, I\ and I 2 can now be written as 


r 2 - 1 ) 

a Fi(l/2,l + r;2 + r;2K/r e ) 


( 12 ) 


€ 


1 + r 


where 2 Fi is a hypergeometric function which can be easily evaluated using a summation 
term as described in Ref. 9. In order to determine the electric field from Eq. (9), the two 
unknowns, and r, must be evaluated. 


C. Ion Current 

In this section we derive a relaxation equation to determine the ion current oscillations. 
We begin with the momentum equation for the ions which also can be written as 


dnu dnu du e dV 

-r— + u— h nu— = n— 

dt dx ox m ox 


(13) 


On multiplying the above equation by unit charge and using the differential of Eq. (5), we 
get 


dl 

dt 


+ 


dJi 

u-r — (- enu 
dx 


du 

dx 


dV 
m, n dx 


( 14 ) 
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In order to obtain an ordinary differential equation in time, we express the spatial deriva- 
tives of u and J t in terms of the electric field. This is done using the derivatives of Eqs. (6) 
and (11). 

du _ _e_dV_ _ _e_V_dV_ 

U dx m dx mV dx 

dx 7 dx v dx 

From Eqs. (14) and (15), we get the following relaxation equation for the ion current, 


dJ± 

dt 


urJi 

~ 


/ V\e 2 idV 

+ V 1 v)m n \dx 


(16) 


We now identify that (u/V)dV/dx /v 1 /r<, where r* is the ion transit time through the 
sheath. This allows us to simplify this equation further. 


where, 


dJj _ (r + g) Jj + 9 ^q^b 

dt ~ Ti 


(17) 


e(V - V) V — V 


9 — 12 


n- 2V 


(18) 


The parameter g is the measure of ion nonequilibrium in the sheath and is the driving force 
responsible for ion flux oscillations. When the ions instantaneously adjust to a relatively 
slow potential variation (i. e. V ~ V^) , g becomes small and ion current oscillations become 
minimal. If g = 0 there is no driving force and the ion current oscillations decay. On 
the other hand when the sheath potential changes rapidly, the high relaxation time, t<, 
prevents the ion flux to respond to the time scale. Thus, Eq. (17) satisfies the general 
properties of ion current relaxation process in a non-collisional sheath. 


D. Solution Technique 

In order to close the system, we also need to determine the value of r. It is emphasized 
that r changes over time, however, its exact dependence on time is unclear. After some 
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numerical experimentation, we choose r — ToV w {t). This assumes that r is proportional to 
the potential, V w , at the electrode over the rf cycle. The constant, r 0 , will be determined 
iteratively from the periodicity of J t , which can be written as 

o (i9) 

Jo dt 

The following steps describe the iterative solution technique once no and T e are obtained 
from the bulk plasma conditions. 

Step 1. Assume initial values of J», V, and ro. 

Step 2. Integrate Eqs. (1) and (17) using a time marching scheme to determine new values 
of Ji and V. ro stays constant until the end of the rf cycle. 

Step 3. Repeat Step 2 until the end of the cycle. 

Step 4. Compute the new value of ro from Eq. (19). If the value of ro is converged, go to 
the next step. Otherwise go to Step 2. 

Step 5. Compute ion velocity and the IED from the values of V, and r over the rf cycle. 
It usually takes less than 25 cycles to obtain converged values of ro- Figure 1 shows that the 
convergence rate of ro is slower at higher values of /3. It is also found that the magnitude 
of the applied potential or current at the electrode does not affect the convergence rate. 

IV. RESULTS AND DISCUSSION 

In this section we demonstrate that the proposed sheath model with ion flux correction 
shows improvements in the prediction of the electric field and ion impact energy distribu- 
tion (IED) at the electrode. We will make comparisons between the present model and 
the unified sheath model. It must be noted that the unified sheath model is obtained by 
simply setting I 2 — 0 in Eq. (9). These models will also be validated against the numerical 
solution of the ion fluid equations, which is considered exact in the current analysis. 

We will consider two cases; potential control, where a time varying potential waveform 
is applied as a boundary condition at the electrode, and current control, where a time 
varying current is passed through the sheath. The results will be analyzed for different 
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values of the frequency ratio (0 = u/u)pi) at T e = 2.0eV with an ion molecular weight of 
69 (CF+). 

A. Potential Control 

We apply a sinusoidal potential waveform at the electrode end of the sheath, 


V w ( t ) = V ac sin u)t + Vdc (20) 

where V dc is the dc self bias which is determined by enforcing a net zero current over 
one rf cycle. The potential is held at zero at the sheath-presheath boundary. A V ac of 
100 V is chosen for all cases run in this section. In Fig. 2 we present the ion current 
oscillation waveform obtained by the present model and the numerical model. A constant 
ion current assumed in the unified sheath model is also plotted. The ion conduction current 
is normalized by the ion current entering the sheath from the plasma. There is reasonable 
agreement between the present model and the numerical solution. The agreement is better 
at 0 = 0.1 and 0.25 than at higher values of 0. The overall trends are well captured and 
the remaining disagreement is attributed to the assumptions made in Eqs. (1) and (11). 
A further refinement of Eq. (17), which is the subject of future study, could remove some 
of the existing discrepancies at higher 0. 

Figure 3 shows the variation of the electric field over the rf cycle. It is evident that the 
ion flux correction improves the prediction of the electric field values. The amplitude 
of the oscillations obtained from the unified sheath model is around 25% smaller than 
the numerical value at 0 — 0.1 and 0.25, while the present model consistently shows 
improved predictions. At 0 = 1.0, all three predictions come closer due to a smaller ion 
flux oscillation. At even higher values of 0 the correction term is expected to be minimal 
and, hence, the use of the unified sheath model will be adequate. 

Figure 4 shows the IED at the electrode obtained from the numerical and the semianalytical 
models. It is important to determine the IEDs accurately since they play an important role 
in the etching process. These distributions are obtained by sampling the ion kinetic energy 
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over the rf cycle into discrete energy bins. The size of the energy bin determines the noise 
that appears in the distribution. The ions acquire a bimodal kinetic energy distribution as 
they travel through a time varying field before they reach the electrode. The two important 
properties of this distribution are the relative peak heights and the energy gap between 
the peaks. The energy gap is determined by the range over which the damped potential, 
V , oscillates, which is predicted fairly accurately by both models. It can be seen in Fig. 4 
that both the numerical solution and the present model predict a higher high energy peak, 
while the unified sheath model predicts peaks of equal heights. The unequal peak heights 
is caused by the ion flux oscillation which weighs in favor of the high kinetic energy ions. 
During the rf cycle when the sheath potential is at its maximum, the flux of ions is also 
high, and hence, more ions hit the electrode with high kinetic energies. The reverse is 
true when the sheath potential drop is at its minim um. It is also interesting to note that 
at J3 = 1.0, the peak heights are not as different as in the case of f3 = 0.25. There are 
two reasons for this. First the ion flux oscillation amplitude is not as high, and second, 
the phase difference between the ion flux waveform and the potential waveform causes a 
weaker constructive interference. 

B. Current Control 

In this case we pass a sinusoidal current through the sheath with no dc component, 


J = J ac sinut (21) 

The sheath potential is obtained as a part of the solution from the following equation, 


J ac sin u)t = enu + J e + eE 


= enu + J e + e-^-zV 


where J e is the electron conduction current. The values of the potential obtained also 
strongly influences the IEDs. In Fig. 5 we show the sheath potential obtained for /? = 0.25 
with the current amplitudes, J ac = 35 and 50A/m 2 . It is evident that the unified sheath 


11 



model predicts an exceptionally large value of sheath potential drop than what is obtained 
from the direct numerical solution of the ion fluid equations. On adding the effect of ion 
flux oscillations, this unrealistically large potential drop is minimized and the result agrees 
well with the one obtained from the full numerical solution. This significant change in the 
potential drop is due to the change in sheath capacitance. In the present model the sheath 
is allowed to store time varying ion charge unlike the unified sheath model where it is held 
constant due to a constant ion conduction current assumption. The electric field waveform 
also shows similar trends. In Fig. 6 we note that the present model predicts values of the 
electric field that agree with the numerical results extremely well. 

Figure 7 shows the IEDs obtained from the two models and is compared with the direct 
numerical solution results. First we note that in the current control case, the low energy 
peak is higher than the high energy peak due to the fact that during the rf cycle the sheath 
potential drop stays near it minimum for a longer time compared to at its maximum [see 
Fig. 5j. Thus, cumulatively a lot of ions hit the electrode at low energies. Although 
the ion flux oscillations favor the high energy peak, it is not enough to compensate for 
the this effect. However, the ion flux oscillations does effect the relative peak heights 
which is evident from a smaller difference in the peak heights in the present model versus 
the unified sheath model. The major difference is, however, seen in the width of the 
distribution. This is clearly an artifact of the large sheath potential drop predicted by the 
unifi ed sheath model, which results is a wider IED spectrum. 

V. CONCLUDING REMARKS 

We propose a semianalytical ion dynamics model that accounts for ion current oscillations 
at intermediate bias frequencies (/? ~ 1). A relaxation equation for ion current is obtained 
which is integrated over the rf cycle along with the equation for the damped potential to 
determine ion flux, kinetic energy, and sheath potential or current. A correction term is 
also proposed for the electric field expression. We see marked improvements in electric 
field and IED predictions when we include the ion current oscillation effect. The biggest 
improvement is found in the prediction of the sheath potential drop when a specified rf 
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current is passed through the sheath. This is due to the change in the sheath capacitance 
since, in this model, the accumulated ion charge is allowed to vary with time in the sheath. 
The relative peak height and the width of the IED are also improved significantly. 
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FIGURE LEGENDS 


FIG. 1. Convergence of tq with rf cycle at different values of frequency ratios. 

FIG. 2. Comparison of normalized ion current, enu/eriQUB , oscillations at V ac = 100 V, 
T e = 2.0eV, and /? = a) 0.1, b) 0.25, c) 0.5, and d) 1.0 A potential control boundary 
condition is used. 

FIG. 3. Comparison of electric field oscillations at Vac = 100F , T c = 2.0eV, and /3 = a) 
0.1, b) 0.25, c) 0.5, and d) 1.0. A potential control boundary condition is used. 

FIG. 4. Ion energy distribution functions at V ac = 100V , T e = 2.0eV , and (3 = a) 0.25 and 
b) 1.0. A potential control boundary condition is used. 

FIG. 5. Comparison of sheath potential at /? = 0.25, T e = 2.0eV, and J a c = a ) 35 and b) 
50 A/m 2 . A current control condition is used. 

FIG. 6. Comparison of electric field at /3 = 0.25, T e = 2.0eV, and J oc = a) 35 and b) 50 
A/m 2 . A current control condition is used. 

FIG. 7. Ion energy distribution at /? = 0.25, T e = 2.0eV, and J oc = a) 35 and b) 50 A/m 2 . 
A current control condition is used. 
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FIG. 2. b) 
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FIG. 2. c) 
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FIG. 2. d) 
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FIG. 3. b) 
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FIG. 3. c) 
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FIG. 3. d) 
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